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Abstract 

The standard model of large scale structure is considered, in which the structure origi- 
nates as a Gaussian adiabatic density perturbation with a nearly scale invariant spectrum. 
The basic theoretical tool of cosmological perturbation theory is described, as well as the 
possible origin of the density perturbation as a vacuum fluctuation during inflation. Then, 
after normalising the spectrum to fit the cosmic microwave background anisotropy measured 
by COBE, some versions of the standard model are compared with a variety of data coming 
from observations of galaxies and galaxy clusters. The recent COBE analysis of Gorski and 
collaborators is used, which gives a significantly higher normalization than earlier ones. The 
comparison with galaxy and cluster data is done using linear theory, supplemented by the 
Press-Schechter formula when discussing object abundances of rich clusters and of damped 
Lyman alpha systems. By focussing on the smoothed density contrast as a function of scale, 
the observational data can be conveniently illustrated on a single figure, facilitating easy 
comparison with theory. The spectral index is constrained to 0.6 < n < 1.1, and in parti- 
cle physics motivated models that predict significant gravitational waves the lower limit is 
tightened to 0.8. 

1 Introduction 

By 'large scale structure' one means galaxies and clusters, with emphasis on their spatial 
distribution and motions, and also the cosmic microwave background (cmb) anisotropy. 
We are at present on the verge of a quantum leap in our understanding of large scale 
structure, because the cmb anisotropy is being measured for the first time |l|, ||, ||, Q. In 
particular, we may in the forseeable future verify or rule out the standard model of structure 
formation, according to which large scale structure arises from a Gaussian, adiabatic density 
perturbation that is nearly scale invariant at horizon entry. If that model is verified there will 
be the dazzling prospect of a window on the fundamental interactions on scales approaching 



the Planck scale, because the vacuum fluctuation during inflation generates just such a 
perturbation. 

This article has two objectives. One is to equip the non-specialist with a starting 
point, from which to follow the saga that will unfold during the coming years. To this 
end we include an extensive discussion of linear perturbation theory, which allows one to 
translate theoretical input coming from say a model of inflation into a form amenable 
for comparison with observations. Our other aim is to provide a critical assessment of 
the present position, in the light of all relevant observations, including the abundance of 
damped Lyman alpha systems and the latest analysis of the crucial COBE data. These 
and the other relevant observations are summarized on a single plot. On the same plot are 
drawn the canonical version of the standard model, and variants which alter the Hubble 
constant h, the spectral index n, and the fraction Q u of any hot dark matter. We discuss 
the observational constraints on these three parameters, and in particular the constraint on 
n, which is of great interest because it can be a sensitive discriminator between inflationary 
models. As we have emphasised several times in earlier publications [||, ^, [7|, ||, ||, 10] a 



powerful constraint on n is provided by the long 'lever arm' between the very large scale 
explored by the large angle cmb anisotropy and the smaller scales explored by galaxy and 
cluster data. On the basis of the presently available data we find that 0.6 < n < 1.1, with 
lower limit tightened to 0.8 in particle physics motivated models of inflation that generate 
significant gravitational waves. The lower limit in particular is rather firm because several 
different types of observation confirm it. 

Three possible models for large scale structure 

Three possible models of large scale structure are commonly entertained. The standard 
model is that it originates as a Gaussian adiabatic density perturbation, whose spectrum 
is nearly scale independent at horizon entry. This model has been explored far more thor- 
oughly than the other two, and is the only one that will be considered here. An alternative 
is that large scale structure originates from topological defects, such as cosmic strings. In 
both of these models the underlying scale invariance means that galaxy and cluster forma- 
tion is directly related to the magnitude of the large scale cmb anisotropy. It also means 
that the first objects to form are the progenitors of galaxies. Finally there is the possibility 
that large scale structure originates from a density perturbation (either adiabatic or isocur- 
vature) with a spectrum that is not even approximately scale independent. In that case 
structure formation could be very different, with perhaps much lighter objects forming first, 
and it would not be directly related to the magnitude of the large scale cmb anisotropy. 

The standard model has two features which distinguish it from the others, and make it 
so attractive. One of them concerns theory. If, as is widely supposed, the initial conditions 
for the hot big bang are set by inflation, then an adiabatic, Gaussian, more or less scale 
invariant density perturbation is predicted. What could be more natural than to suppose 
that perturbation will explain large scale structure? If that turns out to be so, large scale 
structure will provide us with a unique window on the nature of the fundamental interac- 
tions, because both the magnitude and precise scale dependence of the perturbation are 
highly model dependent @, [ll|, IJ, g. 



The other feature concerns phenomenology. The model has been intensively studied 
and is relatively simple, so that by now one knows how to estimate its predictions for most 
available types of data. To a first approximation the predictions depend on a single number, 
specifying the magnitude of the density perturbation at horizon entry. If that number is 
chosen to fit the cmb anisotropy, all other data can certainly be explained to within a 
factor of two or three! Of course the burning question is whether the data can actually be 
explained within their observational uncertainties, which in the best cases are only tens of 
percent. The answer to that question depends what other parameters are available in the 



standard model, and it will be our main focus. 



An overview of the standard model 

The standard model assumes that large scale structure originates as an adiabatic, Gaussian 
density perturbation whose spectrum is more or less scale-independent at horizon entry. 
The scale dependence (if any) is parameterised as a power law, with a spectral index that 
by convention is defined so that n = 1 corresponds to scale invariance. 

In order to have any chance of agreeing with observation, the standard model requires 
non-baryonic dark matter, which is more or less cold, and which has a density dominating 
the baryon density. By 'cold' one means that the constituent particles are stable, non- 
interacting and non-relativistic, at all relevant epochs. 

The simplest version of the standard model assumes that n = 1, and that the (non- 
baryonic) dark matter is completely cold. It also assumes that the energy density of the 
universe is critical, = 1, with no cosmological constant or other exotic contribution so 
that there is the standard cosmology with matter domination at present. This critical 
density, n = 1 CDM model was the favoured one for many years. It contains only two 
free parameters, which are the normalisation of the spectrum, and the value of the Hubble 
constant Hq = 100/ikmsec -1 Mpc -1 . According to observations having nothing to do 
with large scale structure (namely, direct observations and measurements of the age of the 
universe), 0.4 < h < 0.6.[] Assuming the central value h = 0.5 we arrive at a canonical 
version of the CDM model, whose only free parameter is the magnitude of the scale invariant 
density perturbation. 

The canonical CDM model does not agree with observation because the predicted spec- 
trum of the density perturbation has the wrong scale dependence. On large scales the 
spectrum is accurately determined by the COBE measurement of the cmb anisotropy, and 
as one goes down in scale through the regime explored by data on galaxies and clusters it 
becomes progressively too big compared with the data. How can we reduce the small scale 
power? 

The simplest possibility is to reduce h below the canonical value h = 0.5, which reduces 
the small scale density perturbation because it delays matter domination giving the density 
perturbation less time to grow. It has recently been noted jnj| ] that this 'old universe' option 
might work if h is as low as 0.3, but such a value is difficult to reconcile with measurements 
of h from Hubble's law. 

Another possibility is to reduce n below the canonical value n = 1. This 'tilted spectrum' 
option has been widely investigated 0,0, 14, [H], [L6|, [?], |j| and it might be viable with 



n — 0.7. Unfortunately most inflation models with tilt also generate significant gravitational 
waves which ruin this concordance, as we discuss later. Of course one can combine the old 
universe and tilted spectrum options. 

A third option is to change the hypothesis of completely cold dark matter, which reduces 
the small scale density perturbation because cold dark matter maximises its rate of growth. 
An attractive possibility, from both a theoretical and observational viewpoint, is to invoke 
a fraction Q, v of hot dark matter in the form of massive neutrinos. This mixed dark matter 
(MDM) model has been widely investigated f7|, ||, ||, fl], ||, g ||, g(J, both with 



and without the option of allowing h and n to depart from their canonical values. We shall 
see that it may be observationally viable with Vt u = 0.15 or so and the canonical h and 
n. Other possibilities yet to be investigated fully are to replace the CDM by some form of 
warm dark matter like sterile neutrinos, or to replace it by decaying or self-interacting dark 
matter. 



1 Higher values of h are permitted only in low density models, especially those featuring a cosmological 
constant. Given h, the baryon density is practically fixed by the standard nucleosynthesis relation fish 2 = 
0.013 ± 0.002. 



The final possibility is to reduce the matter density, either with [^J, |25|, 26 1 or without 
[27, 25, 26] a cosmological constant to keep the total energy density critical. In practice one 
takes the non-baryonic dark matter to be completely cold in that case. This low density 
CDM model has been quite widely investigated and with h and n at their canonical values 
it may be observationally viable with f2 c of order 0.5 or so. 

All of this assumes that there is no significant gravitational wave contribution to the cmb 
anisotropy. A contribution up to 50% or so is not ruled out by present data, and is actually 



predicted by some models of inflation |j, 11, 12]. Of the models that are well motivated 



from particle physics, those predicting a significant contribution also predict a spectral index 
n < 1, and on the large scales explored by COBE they predict that the relative gravitational 
wave contribution to the mean square cmb anisotropy is R ~ 6(1 — n). In these models the 
normalization of the density perturbation is therefore reduced by a factor [1 + 6(1 — n)] -1 / 2 , 
and the possible existence of this factor should be taken into account when considering tilt 
in the spectral index. On the other hand, almost all inflation models suggest that if the 
spectral index is very close to 1 then the gravitational wave contribution will be negligible. 

In this article the focus is on the (critical density) MDM model, which of course includes 
the critical density CDM model as a special case, allowing h and n to vary, and keeping in 
mind the possibility of gravitational waves. The low density CDM model will be mentioned 
only briefly, and other possibilities not at all. 



2 Cosmological perturbations 

In this section we give a theoretical overview of cosmological perturbation theory, which will 
later be amplified and compared with observation. We also derive briefly the inflationary 
predictions for large scale structure, which will provide a unique window on the fundamental 
interactions if the standard model is verified. 

2.1 Linear cosmological perturbation theory 

The foundation for the standard model of structure formation is linear cosmological pertur- 
bation theory, according to which the perturbations satisfy a set of linear partial differential 
equations as long as they are sufficiently small [28, |29|, |(], [H], 32, 33, 34, [35|, ^g, |37|, 38]. 



Using linear theory, one can follow the growth of an enhancement in the mass density until 
it begins to collapse to form a gravitationally bound structure. One can also give an almost 
complete description of the cmb anisotropy. 

The perturbations decouple into three separate modes, traditionally called scalar, vector 
and tensor modes. Vector perturbations are associated with the vorticity of the fluid flow, 
and are presumably negligible since the vorticity decays. Tensor perturbations are, for 
practical purposes, freely propagating gravitational waves, and so are very easy to handle. 
They are predicted at some level by inflation, and might affect the cmb anisotropy. 

This leaves the scalar perturbations. They are completely defined by giving, for each 
particle species, the perturbation in the momentum distribution function) 2 ]. From this one 
can calculate the perturbation in the matter density, pressure and anisotropic stress of each 
species. (It turns out these in turn determine the metric perturbation, but the latter is not 
directly observable.) 

For a particle species whose collisions are sufficiently frequent to make it a perfect fluid, 
the full momentum distribution is not needed and one needs only the energy density and 
pressure. This is a useful approximation for baryons, electrons and photons before photon 



2 In principle the perturbations in the momentum distribution functions have vector and tensor modes as 
well as the scalar one being considered here, but the vector mode is supposed to be absent, and the tensor 
mode is negligible because gravitational waves have negligible coupling to matter. 



decoupling. Alternatively, if a particle species has negligible random motion it is automati- 
cally a perfect fluid with negligible pressure, and one then needs only its mass density. This 
is the case for cold dark matter and for baryons well after decoupling. These perfect fluid 
approximations provide roughly correct results for both the matter and the radiation, but 
are not adequate to address the increasingly accurate observations now becoming available. 

In order to define a momentum distribution function one needs a set of worldlines, 
specifying the observers who measure the momentum. Also, to define the perturbation 
in this or any other quantity one needs to slice spacetime into hypersurfaces so that the 
quantity can be split up into an average plus a perturbation.^ Any choice of worldlines and 
slicing will do, but a convenient one that we adopt here is to use comoving observers and 
the hypersurfaces orthogonal to them, which are called comoving hypersurfaces. (Comoving 
observers by definition move with the total energy flow, which means that they measure 
zero total momentum density.) 



Independent Fourier modes 

It is convenient to expand each perturbation as a Fourier series in comoving coordinates, 
because in linear theory each comoving Fourier mode evolves independently of the others. 
Comoving Cartesian coordinates x are related to physical Cartesian coordinates r by r = ax, 
where a(t) is the scale factor of the universe. It will be convenient to normalise a = 1 at 
present. A generic perturbation will be denoted by /, and the Fourier series is defined in a 
comoving box much bigger than the observable universe, 

/(x)=£/ k exp(ik.x). (1) 
k 

For each Fourier mode the physical wavenumber is k/a (where k = |k|), and its present 
value is just k. At any epoch, a feature in a perturbation / with size aR is described by 
modes with wavenumber of order ~ a/k. 



Before horizon entry 

A crucial epoch for each scale is horizon entry, when the inverse wavenumber a/k first 
falls within the Hubble distance H~ . For critical density, the scale entering the horizon 
at matter-radiation equality is k^ 1 = (20h~ 1 )h~ 1 Mpc. Smaller scales enter the horizon 
during radiation domination. 

Before horizon entry there is no time for causal processes to operate, and each comoving 



region of the universe evolves like a separate Friedmann universe [35, 37]. The evolution 
of the perturbations is therefore very simple, and can be calculated just by comparing the 
independent evolution of the separate Friedmann universes. 

From now on, we specialise to the standard model of structure formation. The standard 
model states that the initial perturbations are adiabatic and hence determined by the cur- 
vature perturbation, and it also states that the curvature perturbation is Gaussian and has 
a more or less scale independent spectrum. 



The adiabatic initial condition 

The adiabatic initial condition states that well before horizon entry the momentum distri- 
bution function of each species is a function only of the total energy density, or equivalently 
of the temperature. For a perfect fluid species it is enough to specify the energy density 
and pressure, and for a non-relativistic (matter) species the mass density is enough. What 

3 The only exception is if the unperturbed quantity (ie., the average) is time independent, in which case 
the perturbation is independent of the slicing to first order. 



functions of total energy density these quantities are depends on the assumed cosmology, 
the usual assumption being that radiation (extreme relativistic) species are in thermal equi- 
librium and have practically zero chemical potential. 

The adiabatic initial condition determines all the initial perturbations in terms of the 
energy density perturbation. The perturbations are indeed adiabatic, because the evolution 
of every quantity along each comoving worldline is that of a Friedmann universe which is 
adiabatic. (This evolution defines the perturbations, because they all vanish on a hyper- 
surface of constant energy density.) From the relations p m oc a 3 and p r oc a 4 for radiation 
and matter in a Friedmann universe, one learns that the density contrasts of the radiation 
and matter species are related by 

8p m _ 3 Sp r 
Pm 4 p r 

In a Friedmann universe the curvature (measured in a comoving distance unit) is con- 
stant. One can therefore define a curvature perturbation, which is constant on each scale 



well before horizon entry [35]. It is conveniently specified by the quantity 

n = -HSt (3) 

where 5t is the displacement of the comoving hypersurfaces from flat hypersurfaces @ . On 
each scale IZ^ is constant well before horizon entry, and it determines the density pertur- 
bation 5pk. At horizon entry Sp^/p ~ TZ^, so one can think of the curvature perturbation 
as specifying the magnitude of the density contrast at horizon entry. 



Gaussian perturbations 

In comparing theory with observation, one is interested only in the stochastic properties of 
the perturbations. According to the standard model the initial curvature perturbation is 
Gaussian, which means essentially that its Fourier components have random phases (apart 
from the reality condition /k = /* k ). By virtue of the adiabatic condition, this Gaussian 
property is bequeathed to all of the perturbations. 

All stochastic properties of a Gaussian perturbation are determined by its spectrum, 
defined as the modulus squared of its Fourier component. To be precise, we define the 
spectrum Vf of a generic perturbation / as 

V f (k) = ^(Lk/27rf(\f ll \ 2 ) (4) 

where L is the comoving size of the box used for the Fourier expansion, and () indicates that 
|/k| 2 has been averaged over a small region of k space to make it smooth. The spectrum is 
independent of the direction of k because no direction in space is preferred. The normal- 
ization is chosen so that according to Parseval's theorem the mean square perturbation cr 2 - 

V f {k)dk/k . (5) 

Evaluated at random positions, each perturbation has a Gaussian probability distribution 
with variance o" 2 . The spectrum Vt(1) of the cmb anisotropy is defined in an analogous 
way, as essentially the modulus squared of its Ith multipole. 




Scale independence of the initial spectrum 

The linear evolution gives all of the spectra in terms of the spectrum V-n of the initial cur- 
vature perturbation. According to the standard model Vn is more or less scale dependent, 
any scale dependence being parameterised by a power law V-r, oc A; n_1 . (The appearance 
of n — 1 instead of n is a historical accident.) Exact scale independence corresponds to a 



spectral index n = 1, whereas n < 1 corresponds to a 'tilted' spectrum which has less power 
on small scales. One can also contemplate a 'blue' spectrum |], |(], 41], n > 1, but this 
possibility is not favoured by either theory or observation. 



The transfer functions 

Except on very large scales, the perturbations evolve in a complicated manner after hori- 
zon entry, as the particles move about under the combined effect of gravity and particle 
collisions. By the present epoch however, the situation has again become simple. On cosmo- 
logically interesting scales the scale dependence of the matter density perturbation becomes 
fixed, and so does that of the anisotropy of the cmb. The spectrum of the density contrast 
is related to Vn by a linear transfer function, and so is the spectrum of the cmb anisotropy. 
The transfer functions depend in general on the value of h, on the nature and the amount 
of the non-baryonic dark matter and on the value of the cosmological constant if one is 
introduced. The transfer function for the cmb anisotropy also depends on the epoch of 
re-ionisation if it is late. 

Observable scales 

Observations of large scale structure probe a range of scales. 

The scales probed by the cmb anisotropy are easy to estimate. It originates at the 
surface of last scattering, whose distance is now close to 2HQ 1 = 6000/j- 1 Mpc, and the 
thickness of the last scattering surface is about 7/i _1 Mpc. This means that in round figure 
it probes scales 10 Mpc < k^ 1 < 10 4 Mpc. 

Working out the scales that are probed by galaxy and cluster observations is more 
difficult. An upper limit of order 100 Mpc comes just from the fact that much more distant 
regions of the universe have yet to be surveyed in detail, but to describe the lower limit we 
have to describe structure formation. 

In viable versions of the standard model, copious structure formation begins only at a 
redshift of a few, when objects with a very broad range of masses form at more or less the 
same epoch. In order for an object to be reasonably stable, the baryons within it have to be 
able to collapse and dissipate their energy, which is thought to be possible only for masses 
bigger than 1O 6 M or so. The resulting objects are thought to become early galaxies. 
Afterwards there is a bottom-up picture of structure formation, successively more massive 
objects forming in sequence ending at the present epoch with rich galaxy clusters of mass 
around 10 15 Mq. The lower limit of the mass range to which the bottom-up picture applies 
is not very well defined but something like lO n M is a reasonable order of magnitude to 
have in mind 

In the bottom-up picture, each object originates as a slightly overdense region in the 
early universe, the region attracting more matter until it eventually collapses under its own 
weight. As long as its over density is small, the region expands with the universe and is 
described by linear theory. The Fourier modes describing it have inverse wavenumber a/k 
of order its size and in a critical density universe the mass contained in a sphere with radius 
a/k is 

M(k- 1 ) = 1.16 x lO 12 /i 2 (A;" 1 / 1M Pc) 3 M . (6) 

Thus, observations relating to the bottom-up picture of structure formation explore scales 
in the range IMpc < k^ 1 < 10 Mpc. 

What are the relevant observations? Since big clusters formed recently (in the standard 
model), each of them can probably be identified with one of the first objects that formed 
in this mass range. The same may be true of smaller clusters and 'groups', with masses 
1O 13 M < M < 1O 14 M . What about smaller masses? Present day galaxies have 1O 6 M < 
M < 1O 12 M , but because of merging and other astrophysics they can probably not be 



identified with the first objects forming in this mass range. However one also observes 
quasars and damped Lyman alpha systems out to a redshift of 3 or 4. Their masses are 
thought to be in roughly the range M@ ~ W n M & to 1O 13 M0, and they can probably 
be identified with the first objects forming in this mass range. The formation of the first 
objects with mass below this range is not directly observed, but observational constraints 
on the epoch of re-ionisation should provide an indirect probe in the forseeable future. 

The conclusion is that one can observe the formation of the first structure with mass 
from 10 15 Mq down to perhaps 10 11 M Q , thereby exploring the density perturbation on scales 
IMpc < k^ 1 < lOMpc. Bigger scales than this can be explored by looking at the spatial 
distribution and motion of galaxies and clusters. Smaller scales cannot be explored through 
present observations because one does not understand galaxy evolution, but in the near 
future the cmb anisotropy may constrain the epoch of re-ionisation and provide an indirect 
probe on the scale .01 Mpc. 



2.2 The inflationary prediction for the initial spectrum 

The original motivation [^] for scale independence was simply that it is rather natural, 
particularly in view of the fact that V-jzik) is essentially the spectrum of the density contrast 
at the epoch of horizon entry, when non-trivial evolution begins. If inflation sets the initial 



conditions, approximate scale dependence is however predicted [J43| , 35, 44, 45, f46[ . Let us 
see briefly how this works. 

During inflation the energy density is supposed to be dominated by the scalar field 
potential V((f>), and Q is driven to 1 so thatQ 
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where m 2 pl = G^ 1 is the Planck mass and as usual we are setting c = h = 1. The spectral 
index depends only on the first and second derivatives of V which are conveniently defined 
by the dimensionless parameters 



' s S [v) < 8 > 
m pl v " f<i\ 

where primes denote derivatives with respect to 4>. Inflation typically occurs only if 

e<l, H<1. (10) 

This implies that H is slowly varying (on the Hubble timescale) so that 

a oc exp(Ht) . (11) 
Independently of its initial value, <fi typically settles down to the 'slow-roll' approximation 

3H<j> = -V' . (12) 

On any spatial hypersurface the gradient of 5(f) gives the momentum density, so 8<p 
vanishes on the comoving hyper surf aces. The curvature of these hypersurfaces is defined in 
terms of their displacement 8t(x) from flat hypersurfaces by Eq. (g). On flat hypersurfaces 
the inflaton field perturbation is therefore 8(f) = —<f>5t = (fflZ/H. 

In typical models of inflation Sep has negligible interaction with itself and other fields, so 
that each Fourier mode evolves independently (consistent with the assumption that linear 

4 Einstein gravity is assumed, which can usually be achieved by a suitable definition of the relevant fields. 



V 



cosmological perturbation theory applies during inflation). Quantization is straightforward 
and the mass of 8(f) is typically negligible. In order for inflation to be useful cosmologically 
interesting scales have to start well within the horizon, and the corresponding Fourier modes 
of every field must be in the vacuum or the corresponding particles would dominate the 
energy density and spoil inflation j7]]|] The vacuum expectation value of |5</>k| 2 is there- 
fore given by the usual flat spacetime calculation (the scalar field version of the Casimir 
effect that is experimentally verified for the electromagnetic field). It can be evolved for- 
ward in time using the classical equation of motion that 8(f>k satisfies in the Heisenberg 
representation, and using Eqs. (fiol) and (|l^) one obtains the result [44, gg, 45, 46| 



~ H/2n . (13) 

Since H and <fi are slowly varying it follows that the curvature spectrum a few Hubble times 
after horizon exit is given by 

where the star denotes the epoch of horizon exit aH = k. This value is retained until the 
scale k enters the horizon again, after inflation has ended. Using Eq. (|7|) it can be written 

vm^^(^\ v -- as) 
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The observed normalization V^ 2 — 5.8 x 10 5 is a powerful constraint on models of inflation. 

In particular, since e* <C 1 the energy scale must be less than about 10 16 GeV and 
so must the temperature after inflation. This last requirement makes life quite difficult for 
cosmic string theories of large scale structure though the situation is not hopeless pCfl . 
Using Eqs. (|l0|) and (jl2|) the spectral index n = dlogT^/dlog k is given by |fj|, |46| 

n{k) - 1 ~ 2i h - 6e* . (16) 

In typical models of inflation r] and e vary slowly on the Hubble timescale. Since cosmolog- 
ically interesting scales range over only about four decades of k, this means that n(k) can 
be regarded as a constant leading to the power law Viz oc /c n_1 . 

Contrary to the view that is sometimes expressed, there is nothing very controversial 
about this calculation. Provided that 5<j) has negligible interactions its vacuum fluctuation 
well before horizon entry is a standard result of flat spacetime field theory, and its subsequent 
evolution is given by the classical equation of motion. If Eqs. (Ill]) and (12) are valid the 



quoted result is obtained, but if these conditions fail one can still calculate the evolution in 
a straightforward manner using for instance the formalism described in [46|. The predicted 



spectrum will then generally have a strong departure from scale invariance and the model 
will be in danger of being ruled out by observation. Even the assumption of negligible 



interaction can perhaps also be relaxed, leading to a non-Gaussian perturbation [47]. Again, 
the non-Gaussianity endangers the model. 

A similar calculation [^, fl(| for the gravitational waves gives their relative contri- 
bution to the spectrum of the cmb anisotropy, which on large scales is R ~ 12e. For the 
particle physics motivated models in which R is significant, V((f>) is well approximated by 
an exponential or high power-law, which gives r] ~ 2e and hence R ~ 6(1 — n). 

5 We are here making the usual assumption that the curvature is negligible on the present Hubble scale, 
ie., that the density is critical counting the contribution of any cosmological constant. The assumption can 
be relaxed, but then one cannot invoke flat spacetime field theory to define the vacuum. There is however a 
natural definition, and if it is used and the inflaton field is assumed to be in the vacuum the formulas below 
still hold M. 



Measured values of 1 — n or R significantly different from zero would provide a unique 
window on the inflationary potential, and hence on the fundamental interactions. Even the 
upper bounds on these quantities that already exist rule out the usual versions of 'extended 
inflation', because Q such models give both n < 0.7 and R ~ 6(1 — n). 

Now we study in turn the three observable perturbations, which are the matter density 
perturbation, the peculiar velocity and the cmb anisotropy. 



3 The Matter Density Perturbation 

In the Newtonian regime cosmological perturbation theory is a straightforward application 
of fluid flow equations, taking into account where necessary particle diffusion and free- 
streaming (collisionless particle movement). On large scales and in the early universe the 
Newtonian treatment is inadequate and must be replaced by one based on general relativity, 
but here too it is possible to derive fluid flow equations using a treatment strongly resembling 



the Newtonian one fl33, pj, pq, 37 , pq, M . We employ that treatment here instead of the 



more usual metric perturbation approach [28, 31, 32].[] 
3.1 The fluid flow equations 

The strategy is to populate the universe with comoving observers, who define physical quan- 
tities in their own region. A crucial concept is the velocity gradient Uij. It is defined by 
each comoving observer, using locally inertial coordinates in which the observer is instan- 
taneously at rest, as the gradient of the velocity u % of nearby comoving observer, 

= djUi . (17) 

(Note that Ui is defined locally, with only its gradient having global significance.) The veloc- 
ity gradient can be uniquely decomposed into an antisymmetric vorticity ujij, a symmetric 
traceless shear o~ij, and a locally defined Hubble parameter H, 

= H5ij + Oij + Uij . (18) 

In the limit of homogeneity and isotropy, o~ij = = 0. Then the universe is ho- 
mogeneous on comoving hypersurfaces (those orthogonal to comoving worldlines).[] The 
continuity equation is 

^ = -3H(p + p). (19) 

where p is the energy density, p is the pressure and r is the proper time measured by 
comoving observers. The deceleration of gravity is given by 

f = -H>-?f {p + 3p ). (20) 



These lead to the Friedmann equation 



8?rG K 



H< = —?--,. (21) 

The constant K is a measure of the curvature of comoving hypersurfaces, and we assume 
critical density Vt = 1 which corresponds to K = 0. 



The same equations can be derived using the 'gauge invariant' version of the metric perturbation ap- 
proach ]3l], but then the physical significance of the quantities is obscured, and one loses the connection 
with the Newtonian treatment p6[ . 

7 As noted in Section 3 this definition has to be modified if the vorticity is not negligible. 



Now consider perturbations. One can show from angular momentum conservation that 
u>ij decays like (p + p)a -5 . It is therefore presumably negligible though we shall see in the 
next section how it may be handled if desired. We will also see there how to calculate Oij 
which is not negligible, but for the moment we need only H. 

On a given comoving hypersurface, each quantity can be split into an average plus a 
perturbation, 

p(x,t) = p(t)+6p(x,t) (22) 
p(x,i) = p(t)+5p(x,t) (23) 
H(x,t) = B(t)+6H(x,t) . (24) 

The time coordinate t labels the hypersurfaces, and will be taken to be the average of r. 
We would like to choose the space coordinates x = (x l ,x 2 ,x 3 ) to be comoving coordinates, 
related to Cartesian coordinates by r % = ax 1 , with a the average scale factor given by 
a/a = H. This cannot be done exactly, because the expansion is not isotropic and the 
comoving hypersurfaces are not flat. However these effects are of first order, and can 
therefore be ignored when describing perturbations which are themselves of first order. In 
other words all perturbations 'live ' in the unperturbed Friedmann universe. 



Along each comoving worldline the continuity equation Eq. (19) is unchanged. Ignoring 



anisotropic stress the Friedmann equation receives extra terms to become [37 



A JL = + 3p) - ™L + (anisette stress term) . (25) 

dr 3 3 p + p 

This equation is called the Raychaudhuri equation. The anisotropic stress term is typically 
of the same order of magnitude as the exhibited one, or smaller. 

3.2 Separate Friedmann universes 

If the additional terms in Eq. (|25| ) are negligible, the evolution along each comoving worldline 



is identical with that in a Friedmann universe. The Friedmann equation Eq. (21) is obeyed, 
with some K that is constant along each comoving worldline. Assuming critical density its 
average vanishes, so we write K = 5K{x) to remind ourselves that it is a perturbation. 

To see when the additional terms will be negligible it is enough to compare the ex- 
hibited term with the density perturbation — 4ttG5p/3 coming from the second term (it 
is straightforward to verify that the latter does not cancel with other terms, and we have 
just noted that the anisotropic stress term is not expected to dominate the exhibited one). 
For a given Fourier component V 2 = — (k/a) 2 , and since the density is critical the relative 
contribution of the exhibited term is ~ (k/aH) 2 (5p] i /dp] i ). Well before horizon entry the 
adiabatic condition ensures that 5p k < 5pk, and the relative contribution is small. It is also 
small sufficiently long after matter domination, because then 5p becomes negligible. The 
conclusion is that on each scale a comoving region evolves like a Friedman universe before 
horizon entry, and begins to do so again sufficiently long after horizon entry. 

Given this local Friedmann evolution one can calculate the density perturbation in terms 
of the curvature perturbation 1Z ]35| ], It can be shown pl| , [?J that the definition Eq. (|3|) is 
equivalent tq^] 

Perturbing the Friedmann equation gives, to first order, 

8ttG . 5K k 

^-<fyk ; 

6 a 



2H5H k = —5p k - —± . (27) 



8 The factor k 2 means that 1Z k — {?>/2){8K] i /a 2 ){a 2 /k 2 ) measures the curvature perturbation in units of 
the relevant scale a/k. (The intrinsic curvature scalar of the comoving hypersurfaces is R 3 ' = 6K/a 2 .) Yet 
another interpretation of 1Z is that it is essentially the Newtonian gravitational potential caused by 8p. 



We also need to perturb the continuity equation Eq. ([jl]). In doing so one has to remember 
that the comoving worldlines are not geodesies because of the pressure gradient. As a result, 
the proper time interval dr between a pair of comoving hypersurfaces is position dependent. 
One can show (using essentially the Lorentz transformation between nearby observers) that 
its variation with position is given by ]38|, [?| 

dT (28) 



dt \ p + p 

Taking this into account, perturbing the continuity equation gives 

^ = -3(p + p)5H v - 3H5 P]l . (29) 
dt 

Eliminating 5H^ from Eqs. ( |27| ) and (p9|) then gives for the density contrast 5 = 5p/p 
2H~ l d 



5 + 3w dt 



2 



, aH\ 2 2 + 2w 
+ ( "J 4 = (30) 



where w = p/p. During any era when w is constant, Eq. (|30|) has the solution (dropping a 
decaying mode) 



aH\ 2 r 2 + 2w 



kJ s " = 5nto Rk • (31) 



In the radiation dominated era before horizon entry this becomes 

aH\ 2 4 



A; 



J 5 k = g^k (initial) (32) 



and in the matter dominated era it becomes 

(^) 2 <5k = ^ k (final). (33) 

As the labels imply, we are regarding the value of TZ^ during the first era as an 'initial 
condition', which determines its value during the 'final' matter dominated era. 
For future reference note that during matter domination, H oc t~ l oc a~ 3//2 and 

<5k oc a (matter domination) . (34) 

The Newtonian picture 



Eqs. ( |2lD and (|25|), which determine the evolution of the energy density perturbation are 
valid in the framework of general relativity. For the matter dominated case p <C p they have 
however the same form as the Newtonian equations, and at a given epoch the Newtonian 
picture indeed applies to the universe on scales that are sufficiently far inside the horizon [29]. 



Thus we have demonstrated that the Newtonian equations for the density perturbation are 
valid during matter domination, even on scales that are too large for the Newtonian picture 
itself to be valid. We shall shortly demonstrate that the same is true for the equations that 
describe the peculiar velocity. 

3.3 The transfer function 

On scales k^ 1 (20/i _1 )/i _1 Mpc, horizon entry is long after matter domination so that 
the initial and final eras overlap and IZ^ (initial) = 7£k( nna I) • On smaller scales there is a 
transfer function T(k), which may be defined by 



TZ k (final) = T(k)K k (initial) 



(35) 



The corresponding density contrast is given by Eq. (|33|), and it applies to each type of 
matter [p8| . It is Gaussian, and its spectrum is 



p s(k) = ^Vn(k)T\k) . (36) 

The transfer function depends on the nature of the dark matter, and on the value of the 
Hubble constant. The main physical effects are the following [29, 3C, 32] 



• The cold dark matter density contrast grows, as overdense regions attract more matter 
towards them. The growth is slow before matter domination, but rapid thereafter. 

• The baryons, electrons and photons form a tightly coupled fluid before decoupling 
at redshift z ~ 1000, whose density contrast oscillates as a standing wave. The 
oscillations are damped because the photons diffuse out of overdense regions, and 
they carry some of the baryons with them. 

• After decoupling the photons travel freely. 

• After decoupling, the baryon density contrast grows to match that of the cold dark 
matter on scales > 1O 6 M , but continues to oscillate for a long time on smaller scales 
(this may account for the absence of galaxies with M < 10 6 Mq). 

• Massless neutrinos free stream (travel freely) out of any density enhancement. 

• Any hot dark matter (massive neutrinos) free streams while it is relativistic, and then 
its density contrast grows. 

On each scale, the transfer function is applicable after the random particle motion has 
become negligible. In the absence of hot dark matter this happens soon after decoupling 
(on scales bigger than 10 6 Mq). In the presence of hot dark matter it happens only after a 
redshift given by [II]] 

k^ 1 > 0.11(1 + z) l l 2 £l u Mpc . (37) 

As we are interested only in scales > 1 Mpc, the random motion has become negligible by 
the present. It is still marginally significant at z ~ 4 on the scales relevant for the formation 
of the quasars and damped Lyman alpha systems, but the effect is small compared with 
observational uncertainties and we shall not include it. 

The smoothed density contrast 

The linear theory described so far applies only to small perturbations. Except perhaps at 
early times, this means that one has to smooth all quantities on a suitably large scale before 
using it. Adopting the simple 'top hat' prescription, the density contrast <5(x) at each point 
is replaced by its average <5(x, R) within a comoving sphere, whose present radius R defines 
the smoothing scale. Instead of specifying R one can specify the average mass M{R) within 
the sphere, which assuming critical matter density f2o = 1 is given by 

M(R) = 1.16 x W 12 h 2 (R/l Mpc) 3 M . (38) 

Top hat smoothing multiplies each Fourier coefficient by the window function 

It therefore filters out Fourier modes with k^ 1 < R. On the other hand, the spectrum Vs(k) 
of the unsmoothed density contrast decreases rapidly as k" 1 increases on scales k^ 1 > 1 Mpc, 



corresponding to the fact that the probability of finding a feature with size a/k decreases 
rapidly. As a result, most of the structure in the smoothed density contrast at a given epoch 
has size of order aR, for smoothing scales R > 1 Mpc. 

The mean square <J 2 (R) of the smoothed density contrast is 



At any epoch it can can be made arbitrarily small by choosing a sufficiently large filtering 
scale R, but it grows with time. As long as cr(R) < 1, linear theory applies except in those 
rare regions of space where <5(x, R) > 1. These regions correspond to gravitationally bound 
objects with mass of order M(R) and linear theory does not apply to them, but since they 
are rare they will not significantly affect the evolution of cr(R), or of the Fourier modes with 
k^ 1 > R that dominate it. When o~(R) grows to become of order 1 on the other hand, the 
formerly rare gravitationally bound regions become common, and the evolution becomes 
non-linear on scales k^ 1 ~ R. Since o~(R) oc a oc (1 + z)~ l (Eq. (|34|)), this occurs at a 
redshift given by 



where the denotes the present value of the linearly evolved quantity. As we shall see z n \ 
is no more than a few in observationally viable versions of the standard model, even on 
the smallest relevant scales. Before this epoch gravitationally bound objects are rare, and 
whether they are too rare to account for what is observed is one of the crucial issues that 
we shall address. 

4 The peculiar velocity 

Associated with the density perturbation is a peculiar velocity field, specifying the departure 
from uniform Hubble flow that occurs as matter falls into overdense regions. 

4.1 The Newtonian regime 

The Newtonian picture |^9| applies to the universe after matter domination on scales suffi- 
ciently far inside the horizon, and in particular it applies at present to the sphere of radius 
a few hundred Mpc around us within which detailed observations of galaxies and clusters 
have been made . 

In the Newtonian picture there is a globally defined fluid velocity field u, and choosing 
the reference frame so that u vanishes at the origin the peculiar velocity field v is definecQ 




(40) 



1 + z nl (R) = a {R) 



(41) 



by 




(42) 



An equivalent statement in terms of the velocity gradient Eq. (18) is 



Suij = diVj . 



(43) 



where di = a~ l d/dx l . Like any vector field, v can be written v = v L + v T , where the 
transverse part v T satisfies divf = and the longitudinal part is of the form v L = V"i/V 
Eq. ( |l8| ) defines the local Hubble parameter H, shear and vorticity Uij, this last being 
given by 




Angular momentum conservation gives uiij — > a 2 , so v T decays like a 
dropped (remember that di = a^d/dx 1 ). 



and may be 



9 Up to a constant, which can be chosen so that the average of v vanishes. 



Taking v to be purely longitudinal, it is determined by the density perturbation in the 
following way. First, take the trace of Eq. ( |l8|) to learn that V.v = 35H. From Eqs. (P7|) , 
(p6|) and ( [3l] ) and the matter domination relation Ht = 2/3, it follows that 



V.v = -{4nG5p)t . (45) 

The solution of this equation is 

v = -tVip (46) 

or 

« i (x,t) = -(t/o)^^ (47) 

where 

^(x) = -GoT 2 ( ^4<*V . (48) 
J |x x| 

The factor a~ 2 converts coordinate distances into physical distances. Since tp is related to 
the density perturbation by the Newtonian expression, it is called the peculiar gravitational 
potential. It is independent of t because, from Eq. (|34]), 5p oc a 2 . 

From Eq. ( |33| ) we see that the peculiar gravitational potential is related to the spatial 
curvature perturbation by 

V> = --ft (final) . (49) 
5 



From Eqs. (46) and (Eq) the Fourier components of v, ip and 5 are related by 



v k . ,50, 
* k - -1(f)' (51) 
4.2 Peculiar velocity in general relativity 

The Newtonian picture has to be replaced by general relativity before matter domination, 
and even after matter domination on scales that are not far enough inside the horizon. 



To define peculiar velocity in this case one can proceed as follows 0, 36]. Start with the 



definition Eq. (17) of the velocity gradient, and first assume that the vorticity is negligible 
since angular momentum conservation ensures that it decays like (p + p)a~ 5 . Then one can 
show [|36| that the velocity gradient perturbation is of the form 

6uij = diVj + ^hij (52) 

where v is a globally defined peculiar velocity field as in the Newtonian picture and the 
extra term hij is transverse, djh{j = 0, and traceless, 5 l3 hij = 0. The extra term represents 
the effect of gravitational waves, and at any epoch after matter domination the Newtonian 
picture is valid on scales well inside the horizon provided that it is negligible. It is certainly 
negligible at the present time in the region of the universe around us with radius a few 
hundred Mpc, where detailed galaxy surveys are performed, but it could be significant on 



bigger scales and at earlier times and so contribute to the cmb anisotropy 36]. Its presence 
does not however spoil Eqs. (p5|)-(|5"i~|) because it is transverse and traceless. The conclusion, 
as advertized, is that the Newtonian equations of the last subsection apply to the peculiar 
velocity on all scales after matter domination. We shall use them in the next section to 
derive the Sachs- Wolfe effect. 

One can avoid dropping the vorticity by proceeding as follows p6| . First, a transverse 
velocity v T may be defined by Eq. (H3J). The worldlines with velocity — v T relative to the 



comoving worldlines have no vorticity, and comoving hyper surf aces are defined to be orthog- 
onal to them (there are no hypersurfaces orthogonal to worldlines with nonzero vorticity). 
The velocity gradient Suij receives an extra contribution \{divj — djvj) + ^(diwj + djw] 

where w T 



l + 6(l + f) (?f)' 



r . For scales well inside the horizon w T is negligible 

and the Newtonian result is recovered. 

These results demonstrated, within the fluid flow approach, that the perturbations de- 
couple into three modes. One is associated with gravitational waves, another with the 
vorticity of the velocity flow, and the third with the density perturbation that is our main 
concern. 



5 The cmb anisotropy 

The first detection by COBE of the intrinsic anisotropy of the cmb J|], announced in 1992, 
was surely the most important advance in observational cosmology for a long time. The 
COBE satellite explores large angular scales, corresponding to linear scales of order the 
particle horizon distance 2i7 ~ 1 ~ 10 4 Mpc, and if the standard model is correct it provides 
a direct measurement of the primordial curvature perturbation. More recently other obser- 
vations have detected smaller scale cmb anisotropy, which in due course will probe the finer 
details of the standard model and help to prove or disprove it |||, 

The anisotropy of the cmb is defined by giving the variation of its intensity with direction, 
at fixed wavelength. This variation is usually specified by giving the equivalent variation in 
the temperature of a blackbody distribution. [j Denoting the direction of observation by a 
unit vector e, the anisotropy may be expanded into multipoles 

^ = w.e + £ £ ariT(e). (53) 

1=2 m=-l 

The dipole term w.e is well measured, and is the Doppler shift caused by our velocity 
w relative to the rest frame of the cmb. Unless otherwise stated, AT will denote only the 
intrinsic, non-dipole contribution from now on. 

A feature in the sky of angular size 9 radians is dominated by multipoles of order 
/ ~ 1/9. (This is analogous to the relation R ~ 1/k between the linear size of a feature and 
the wavenumbers that dominate its Fourier expansion.) Translating to degrees we have the 
following relation between I and the angular scale 

(54) 

The cmb originates at high redshift, and therefore on a comoving sphere whose present 
distance is close to the particle horizon at 2Hq = 6000/i~ 4 Mpc. Its anisotropy also 
originates at comparable distances ^] A feature with angular size 9 therefore corresponds 
to a linear scale 

9 9 

R ~ 2# r T 1 = 100/r 1 Mpc-pr . (55) 

1 radian l u 

Equivalently, the Zth multipole of the cmb explores the scale 

R ~ 2H 1 /l = 6000/T 1 Mpc// . (56) 

Present observations have resolution of order 1 degree, so they give information on scales 
100 Mpc < k^ 1 < 10 4 Mpc, and the lower limit will drop to around 10 Mpc when the angular 
resolution is improved. There is no significant anisotropy on smaller scales, because it is 
washed out by the finite thickness ~ 7/i _1 Mpc of the last scattering surface. 

10 In principle the anisotropy depends on the wavelength but the dependence is not significant. 
n Except for the Sunyaev-Zel'dovich effect, which is caused by galaxy clusters. It is significant only on 
arcminute scales, and will be ignored here. 



5.1 The spectrum of the cmb anisotropy 




The value of each multipole af 1 depends on the observer's position. For a random choice 
of position there is some probability distribution, with zero mean, and some variance 
(I ™ Imposition wn i cn is the expected value of {af 1 ] 2 for a randomly placed observer. The 
variance is independent of m because no direction is space is preferred. By analogy with 
Eq. (0), one may define the spectrum Vt(1) of the cmb anisotropy by 

?>,(*) -^^(Kf) . (57) 

47T \ ' position 

The normalization is chosen so that the squared anisotropy averaged over both direction in 
the sky and position of the observer is 

oo 

= EW. (58) 

■ - 9 

position 

The cmb anisotropy is well described by linear theory, so the multipoles at a given 
position depend linearly on the Fourier components of the initial curvature perturbation 1Z. 
According to the standard model 1Z is Gaussian, which means essentially that the phases 
of its Fourier components are random. As a result the spectrum Vt(1) depends linearly 
on the spectrum Vfc(k) of the initial curvature perturbation, so their is a 'photon transfer 
function' Tf(k) such that Vt(1) = Jq° 'rfifyT'ntfydk. Unless there is early re-ionisation, 
which we shall see is unlikely in viable versions of the standard model, the transfer function 
is determined by the amount and nature of the dark matter, through the same physical 
processes that we listed already. The only difference is that, except at the bottom end of 
the relevant range of scales lOMpc < k~ l < 10 4 Mpc, it is not sensitive to the nature of the 
dark matter. 

The spectrum Vt{1) gives ( | af 1 1 2 } position; but we can only measure |a™| 2 at our position. 
The best guess is that it is equal to (laj™! 2 position i but one can also calculate the variance 
of this guess, which is called the cosmic variance. Because 7Z is Gaussian, the real and 
imaginary part of each multipole has an independent Gaussian distribution and as a result 
the cosmic variance of J2m \ a T\ 2 ls om y 2/(2/ + 1) times its expected value. However, 
typical present observations suffer from an associated effect, sample variance [jEJ, which 
arises because they cover only part of the sky preventing measurement of the full set of af 1 . 
In essence, cosmic variance is the part of the sample variance which remains even once the 
full sky is observed, due to the recognition that complete determination of the spectrum 
would involve measurement on an ensemble of microwave skies rather than on the single 
one which is accessible. 

A prediction for the spectrum is shown in Figure 1, and a qualitatively similar result is 
obtained in any viable version of the standard model. 

5.2 The large scale cmb anisotropy (Sachs- Wolfe effect) 

The simplest regime is the flat part of the spectrum corresponding to I < 30 or angular scales 
bigger than a few degrees. If the standard model is correct the anisotropy in this regime 
probes directly the initial curvature perturbation, being caused essentially by the variations 
in the peculiar gravitational potential encountered by the cmb on its journey from the last 
scattering surface. This effect is called the Sachs- Wolfe effect, and we calculate it now for 
the case of critical density and zero cosmological constant following the approach of 0, |36| . 
The generalisation to the other cases [24, |25|, 27 1 is completely straightforward, except for 



the lowest few multipoles where the effect of spatial curvature may become important 



It is convenient to separate the anisotropy into two parts, 



*m_r*rm . (59) 

\ J / initial V / jour 

The first term is the initial anisotropy, measured soon after decoupling by a set of comoving 
observers whose detectors point radially outwards and whose clocks are synchronized (this 
last condition just says that the anisotropy is defined on a comoving hypersurface) . The 
second term is the additional anisotropy acquired on the journey towards us, including the 
Doppler shifts associated with the peculiar velocities of the initial comoving observers and 
of ourselves. (Our own peculiar velocity contributes only to the dipole but for the moment 
we are retaining all multipoles.) 

On the scales in question the adiabatic initial condition Eq. (g) still holds at decoupling. 
Since p r oc T 4 it follows that 

(^Q) ^^(x em ,W) (60) 

V 1 / initial «J 

where 5 is the matter density contrast. Here i em denotes the time of emission (decoupling), 
and x em = 2H ~ 1 e denote the comoving coordinates of the point of emission. 

To calculated the second term, consider a photon passing a succession of comoving 
observers j?], j36fl . Its trajectory is adx/dt = — e and between nearby observers its Doppler 
shift is 

dA da . . 

— = eiBjUijdr = — — + eiejduijdr (61) 

A a 

where the first term is due to the average expansion, and the second is due to the relative 
peculiar velocity of the observers. Integrating this expression gives the redshift of radiation 
received by us, which was emitted from a distant comoving source. The unperturbed result 
is A/A em = 1/flonn and the first order perturbation gives 

/AT(e)\ f Xeia 



\ T 



jour 

The integration is along the photon trajectory 



pern 

/ eiejSuijfot, t)a(t)dx . (62) 
Jo 



dt /tn t\ 
x(t)= - = 3 -° -~ . (63) 



it a V°Q a / 

Neglecting any gravitational wave contribution, 5uij is the gradient of the peculiar velocity, 



Eq. (43), which is given in terms of the peculiar gravitational potential by Eq. (f46[). Using 
Eq. (63) and integrating by parts one finds 



(^P-) = e • [v(0, t ) - v(x em , t em )] + \ [^(x em ) - V(0)] • (64) 

V 1 / jour 



The term e.v(0, to) is the Doppler shift corresponding to our peculiar velocity (in linear 
theory). It contributes only to the dipole and will be discounted from now onQ The term 
— e.v(x em , t em ) is the Doppler shift corresponding to the peculiar velocity of the initial ob- 
server. The remaining term therefore gives the anisotropy acquired on the journey towards 
us if the initial and final anisotropics are defined by observers with zero peculiar velocity in- 
stead of by comoving observers. It involves the gravitational potential ifi and can be thought 
of as being due to gravity, though we derived it by considering a sequence of Doppler shifts. 

12 By manipulating this expression it can be shown that in linear theory the dipole is precisely the Doppler 
shift due to our motion relative to the average motion of everything inside the sphere that emitted the cmb 



The term —ip(0)/3 can be dropped since it does not actually contribute to the anisotropy, 
leaving only the term ^(x em )/3. 

Eq. (m) is exact in linear theory, but we are considering only scales that are far outside 
the horizon when the cmb is emitted. On these scales Eqs. ( |50| ) and (|5l|) show that 

W > |<5k| » |v k | . (65) 

As a result we can drop the initial Doppler shift, and we can also drop the initial anisotropy 
Eq. (|60|). This gives the Sachs- Wolfe formula 

AT(e) 1 , 

—^-L = -V(x em ) . (66) 

It gives the large scale cmb anisotropy in terms of the peculiar gravitational potential, which 
is related to the curvature perturbation by tp = —(3/5)7?.. On the large scales that we are 
considering 1Z retains its initial value. 



The Sachs- Wolfe spectrum 

Inserting the Fourier expansion of TZ and projecting out the multipoles one finds the spec- 
trum [51] 



V T {1) 
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jf(2k/H )V n (k) 



(67) 



where ji is the spherical Bessel function and V-r is the spectrum of the initial curvature 
perturbation. For a power-law spectrum V-ji(k) = (2k/Ho) n ~ 1 Vii(Ho/2) one finds 



V T (l) 



1 21 + 1 
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For n = 1 the square bracket is equal to 1. 



xJLin + 1i r((3-n)/2)r(l + (n-l)/2) 
2 1 j r((4-n)/2)r(/ + (5-n)/2) 



Vn(H /2) 



(68) 



The COBE measurement of the cmb anisotropy 

The COBE observations provide (after considerable analysis) estimates of the multipoles 
a™ in the range 2 < I < 30. In this regime the Sachs- Wolfe formula applies leading to 
Eq. (|68|), and the multipoles are found to be consistent with this expression which is a 
highly non-trivial test of the standard model. From Eq. (^) one can estimate the spectral 
index n and the normalization T , tz(Hq/2). Instead of the latter one can specify any of the 
quantities Vt(J), which define the expected squared multipole. The usual quantity is the 
quadrupole, defined by 

QLs-ps = 2T 2 V T (2) . (69) 

From Eq. (|58|), the factor 2T 2 makes Q 2 ms _p§ the quadrupole contribution to the expected 
mean square of AT. 

The latest analysis of the COBE data [Q] finds that a good fit can be obtained with 
the spectral index anywhere in the range n ~ 0.6 to 1.4. At fixed n the normalization is 
however well determined, and in particular if n = 1 then 

Qrms-PS = (19.9 ± 1.6)^K , (70) 

which corresponds to V^ 2 = 5.8 x 10~ 5 . 

The best fit value of Qrms-PS depends significantly on n, but according to the analysis 
of it becomes practically n independent if the 9th multipole is used instead. This fixes 
the normalization for all n. 



5.3 The smaller scale cmb anisotropy 



Once the normalisation has been determined, the smaller scale anisotropy can be calculated 
in a given version of the standard model. It is practically independent of the nature of the 
dark matter except on very small scales, and barring the possibility of early re-ionisation 
is therefore determined by n, h, the matter density £l m and any cosmological constant 
contribution In the (critical density) MDM model that we are focussing on it therefore 
depends only on n and h. 

As accurate measurements of the smaller scale cmb anisotropy are only beginning to be 
available firm conclusions are not yet possible, but this situation will change in the near 
future. 

As shown in Figure 1, the predicted spectrum exhibits a series of peaks and troughs as 
one goes down in scale, which are usually called 'Doppler peaks'. As an approximation, 
they may be regarded as coming from the standing-wave oscillation of the baryon-photon 
fluid density perturbation, between horizon entry and decoupling. The frequency of the 
oscillations goes up as the scale goes down, so the magnitude of the perturbation at decou- 
pling is an oscillating function of scale. It causes both an initial anisotropy ATi n i t i a i as seen 
by a comoving observer, and an additional anisotropy ATj our corresponding to the peculiar 
velocity of this observer. 

Re-ionisation before z = 50 or so would markedly reduce the height of the peaks, but 
as we discuss at the end of the next section such an early epoch is not predicted by viable 
versions of the standard model. The observations seem to be seeing the peaks at more 
or less the expected height Q in weak confirmation of this prediction, but a quantitative 
analysis is not yet possible. 

6 Theory and observation 

Now we return to the density perturbation, and see what constraints are placed on it 
by observations of galaxies and clusters. We focus on a(R,z), the rms of the linearly 
evolved smoothed density contrast evaluated at an epoch corresponding to redshift z. As 
we noted earlier, apart from small corrections if a hot dark matter component is present it 
is proportional to (1 + z) , so it is enough to consider the present value 



Practically all of the data can be presented in terms of ctq(R), which makes it preferable to 
the more widely used power spectrum Vo(k). (A subscript will always denote the present 
linearly evolved quantity.) 

Some data points and theoretical curves are shown in Figures 2a and 2b, the latter being 
a close up of the former. The curves are all normalized to the COBE cmb anisotropy as 
described earlier. The full line corresponds to the canonical CDM model (n = 1, h = 0.5 
and Q u = 0), using the CDM transfer function of reference and taking £Ib = 0.05 
because of the nucleosynthesis relation Hg/j 2 = .013 ± .002. The other lines show the effect 
of varying one at a time the parameters n, h and Vt u . To vary Q u at fixed h = 0.5 we have 
used the parameterisation of the transfer function in reference plj , and to vary h at £l u = 
we have used the approximation that the CDM transfer function depends only on kh~ 2 . 
(On small scales and in the presence of hot dark matter highly accurate transfer functions 
have yet to be published but these are probably adequate for the purpose at hand.) 

Since the curves are rather similar and vary over a couple of orders of magnitude, the 
difference between them is best exhibited by normalising them all to a single model, which 
we take to be the canonical CDM model. The result is shown in Figure 3. In what follows 
the prediction of the canonical CDM will be denoted by o"o anon (i?). 




(71) 



To the extent that all of our COBE normalized curves come close to crossing at R ~ 
4000/1- 1 Mpc, the COBE data can be regarded as an observational point for <7o(4000/j 1 Mpc) 
and for convenience this is shown in the Figures. Its error bar is significant because it shows 
how much the theoretical curves can be moved up or down, but we emphasise that the nor- 
malization of each curve is determined directly from the cmb anisotropy as described above, 
not by fitting to this point. 

The rest of the data points concern galaxies and clusters, as we now discuss. 



6.1 The distribution of galaxies and galaxy clusters 

The most extensive observations concern the distribution of galaxies and clusters, and pro- 
vide information about o~o(R) on scales 10 Mpc < R < 100 Mpc. They also provide infor- 
mation about the correlation function £,o(R). These quantities are related to the spectrum 
by 

al(R) = / W 2 (kR)V (k)dk/k (72) 
Jo 

Co{R) = / W(kR)V (k)dk/k (73) 
Jo 

where W is the 'top hat' window function, and £o(R) is taken to be the volume averaged 
quantity [|5^] . A subscript always denotes the present linearly evolved quantity. 

The number density contrasts are known fairly well for IRAS galaxies, optical galaxies, 
radio galaxies and Abell galaxy clusters, out to a distance of several hundred Mpc. Within 
the rather large observational errors they are consistent with the biasing hypothesis, 

5i(x) <5 (x) <5r(x) S a (x) 

— — ~ — ~ — ~ — ~ e)(x) (74) 

oi O b R bA 

where 5j etc. are the number density contrasts, 5 is the matter density contrast and bi 
etc. are scale independent bias factors close to unity. The subscripts I, O, R and A refer 
to the four classes of objects mentioned. One can use the data to estimate the bias factors 
and hence the linearly evolved smoothed density contrast, on scales 1 Mpc < R < 100 Mpc. 
(The determination of the overall normalization comes from nonlinear corrections, including 
redshift space to real space corrections.) 

From a number of possibilities, we have chosen the analysis of Peacock and Dodds [52], 



which combines a variety of data sets. These authors find bj : bo '■ b R : = 1 : 1.3 : 1.9 : 4.5 
for the ratios of the bias factors, and bj = 1.0 ± 0.2 for the overall normalization. Knowing 
the bias factors, the analyses of various groups gives estimates of Vo(k), o-q(R) and S,o(R)- 
Peacock and Dodds convert the last two into estimates of Vo(k) using the prescription^ 

MR) = vl /2 {k R ) (75) 
UR) = vl /2 {V2k R ) (76) 



where 

k R 



1 (m 



l/(m+3) /F 

-r < 77 » 



.2 V 2 

and m = (k/'Po)(dT , o/dk) is the effective spectral index evaluated in the canonical CDM 
model. These formulae are obtained by taking m constant, and using the approximation 

W{kR) = exp{-k 2 R 2 /10) (78) 

13 Although this prescription relating cro(-R) to Vo(k) is adequate on scales > 10 Mpc, it is too dependent 
on the shape of the transfer function to be useful on smaller scales. As we shall see, data on such scales 
directly constrain ao{R), which is our main reason for focusing on that quantity rather than on "Po(fc) or 



which is exact for kR <C 1. We have used this same prescription to convert the Peacock- 
Dodds estimates of T > o(k) into estimates of ao(R), and the results are shown as open squares 
in Figures 2 and 3. For clarity a point with particularly large error bars, which lies between 
the two largest scale points shown, has been omitted. 

Most of the uncertainty in this determination of ctq(R) comes from the uncertainty in 
the overall normalization, defined by bj = 1.0 ± 0.2. If all of the raw data used in the 
determination referred directly to real space number densities, do would be proportional to 
bj times the raw data, so its uncertainty would be the same as for bj, namely 20%. This 
is the case for the results of the APM survey, which is a very important part of the total 
input, but as Peacock and Dodds note the rest of the input is based on redshift space number 
densities. In the linear regime the determination of do from raw data on a given type of 
object in redshift space is proportional to [bi(l + .66rbJ 1 + .20r 2 bJ 2 )] _1 , where r = 6//6j\r- 
For clusters r = 4.5 so that do is more or less proportional to bj 1 as for real space data, but 

—5/7 

for the three types of galaxy r = 1.0, 1.3 and 1.9 and do is roughly proportional to bj . 
Thus, galaxy data in redshift space in the linear regime gives an uncertainty in do(-R) of 
order (5/7) x 20% = 14%. To be on the safe side we work with 20%, noting that at least in 
the linear regime the true uncertainty is between 14% and 20%. 

With bj fixed, one obtains comparatively small error bars, which are indicated by the 
inner error bars in Figures 2 and 3. As an aid to clarity, in Figure 3 we have connected 
the ends of these error bars by straight lines to form a box. Within the box the points can 
be moved up and down more or less independently, and the whole box can be moved up 
or down by 20% as indicated by the outer error bars. Both these types of error are the 
result of statistical fitting, and presumably correspond to something like a 1-d level (with 
the caveat that the outside error bars might be reduced somewhat as discussed above). 



6.2 The bulk flow 

The smoothed peculiar velocity field v(i?, x) is called the bulk flow, and on scales above 
10 Mpc it is described by linear theory. It is therefore the gradient of a potential, and can 
in principle be constructed from its radial component j53|. The latter is directly observed 
from the redshift-to-distance ratios of galaxies, in the region around us with radius a few 
hundred Mpc. 

Unlike the density contrast, the bulk flow can reasonably be assumed to be the same 
as that of the underlying matter at least on large scales, which in principle allows one to 
dispense with the biasing hypothesis. In practice one still needs the hypothesis at present if 
the data are to yield really powerful results, which are obtained by comparing the density 
field obtained from the bulk flow with that obtained from galaxy surveys. A recent study 



[54] concludes that at 95% confidence level bj = 0.7 +0 '" 



0.2- 

Although the bulk flow alone does not yet give very powerful results, it is not completely 
useless. The standard way of utilising it is to look at the magnitude v(R, 0) of the bulk flow 
at our position, or in other words at the average peculiar velocity in the sphere around us 
of radius R-f^\ From Eq. ( ]50| ) the expected value of v 2 (R, 0) is 

al(R) = f (^) 2 W 2 (kR)P (k)f . (79) 

Evaluated at random positions, each of the three components of the bulk flow has a Gaussian 
distribution with one third of this variance, so at 1-d level one expects that 



v{R,0) _ 1+3 2% 
a v (R) ~ " 52% 



US • (80) 



14 There are several variants of this procedure, which correspond to using window functions different from 
a top hat. In what follows we are actually using the variant described in «n. 



Observational values of v(R, 0) are available only for R of order tens of Mpc, and noth- 
ing is lost by considering just R = 40/i -1 Mpc. The canonical CDM model gives || 
cr£ anon (40/i -1 Mpc) = 460 km sec -1 , and the observational value is t> o bs(40/i -1 Mpc, 0) = 
(400 db 15%) km sec -1 . This error is swamped by the 'cosmic variance' given by Eq. (p0|), 
and ignoring it we deduce that 



a v (40h 1 Mpc) _ f400/460)+ 92% 
cT™(40/i -1 Mpc) " ^ 4UU / 4b ^-24% 



(400/460)™ (81) 



where the error is just the inverse of that in Eq. 

Compared with Eq. (40) for ao(R), Eq. (|79| ) contains an extra factor (Ho/k) 2 and as a 



result a v (40h 1 Mpc) probes about the same scale k 1 as eg (90/i -1 Mpc). We conclude that 
the above ratio applies to this latter quantity, which leads to the estimate ao(90/i -1 Mpc) 
.044+;° 4 ° shown in Figures 2 and 3. 



6.3 Abundances 

The rest of our observational points invoke the well known Press-Schechter estimate j/j for 
the contribution to the mass density of objects with mass > M, 

n(> "<*>-">-«*fcsfej) (82) 

where 5 C = 1.7 and M(R) is the mass in a comoving sphere of present radius R, given by 
Eq. (|38[). Except for a more or less unmotivated factor 2 the right hand side is just the 
fraction of space occupied by regions of space where the linearly evolved smoothed density 
contrast exceeds S c , so roughly speaking the Press-Schechter estimate simply states that 
such regions are occupied by objects with mass bigger than M. It is only supposed to be 
valid as long as these regions are rare, corresponding to 0(> M, z) <C 1. In this regime 
cr(R, z) -C 1, and linear evolution should be valid. 

The value 5 C = 1.7 is motivated by a spherical collapse model, and N-body simulations 
suggest that the Press-Schechter estimate is roughly correct, if 5 C is within ten percent or 
so of this value g, 



In some cases the observational result is given directly as a value of Cl(> M), but in 
others it is given as a number density n(> M). In the latter case we obtain an observational 
value of fl(> M) by setting p[> M) ~ Mn(> M) which gives 

n(> M) = 3.60h—^— „ n i* M l , • (83) 
v 1 1O 12 M (/i -3 Mpc) -3 ^ ' 

This overestimates f2(> M) and therefore ct(R), but in practice the mass of the observed 
objects is too uncertain for this to matter very much. 

An alternative is to obtain a theoretical prediction for n(> M) from the Press-Schechter 
formula. This more complicated procedure has been used by many authors (starting with 
the originators of the formula) , but because of the uncertainty in the mass there is generally 
no significant advantage in using it. Some relevant comparisons are cited below. 



Galaxy clusters 

We first look at galaxy clusters. The most useful are those of 'richness class R > 1', 
which are the heaviest bound objects. Their present number density is known to be 
about n ~ 8 x 10 -6 /i 3 Mpc -3 and their mass M c i us is roughly 1O 15 M0. The uncertainty 

15 Somewhat smaller values are found with a Gaussian instead of a top-hat smoothing but we are using 
the latter. 



in the mass is the limiting factor in the present context. Let us first set it equal to 
1.2 x 1O 15 M which corresponds (with h = .5) to R = 8/T 1 Mpc. Then Eqs. @ and 
(H) give o- (8/i _1 Mpc) = 0.71. Lowering M c i us by a factor 2 gives do^Ah' 1 Mpc) = 0.64, 
and raising it by a factor 2 gives ao(W.lh~ 1 Mpc) = 0.80. These three values are indi- 
cated in Figures 2a and 2b by a central point and a slanting 'error bar'. Because OQ anon (i?) 
is strongly varying the corresponding uncertainty in ao(R) / aQ anon (R) is quite big, namely 
+43% and —27%. Compared with this uncertainty the range of scales R = 6.4/i _1 to 
Wh- 1 Mpc is negligible (essentially because all the spectra we consider have broadly similar 
shapes), and we have ignored the horizontal error when plotting Figure 3. 

Our conclusion is that if R > 1 clusters have mass M = 10 15 M Q within a factor 2, the 
normalization of the standard model of structure formation is o"o(8/i _1 Mpc) = 0.71^7% = 
0.5 to 1.0. Estimates of the normalization in a similar spirit have been made by previous 
authors. A recent and typical one |)7j quoted a normalization a(8h~ 1 Mpc) ~ 0.57, but did 
not quantify the effect of the uncertainty in the cluster mass. 

One can in principle apply the same technique to clusters observed at high redshift. A 
treatment of z = 0.3 observations J58| which estimates the effect of both mass and number 
density uncertainty gives o"o(8/i -1 Mpc) = 0.6 to 0.9, in excellent agreement with our result. 
High redshift observations of clusters are likely to become increasingly important in the 
future. 



Quasars and damped Lyman alpha systems 

The same technique can be applied to quasars and damped Lyman alpha systems observed 
at redshift 3 to 4. Their masses are thought to be roughly in the range 10 11 to 1O 13 M and on 
these scales the Press-Schechter formula gives Q(> M, z) <C 1, indicating that linear theory 
applies and the formula is valid. An observational value of the quasar fraction O quas (> M, z) 
or the damped Lyman alpha system fraction ^i ya i p h(> M, z) therefore provides a lower limit 
on f2(> M,z), and on o"o(-R).[3 ^ qU as(> M,z) can be estimated from observation out to 
z = 4, in the range roughly lO n M to 1O 13 M0, if one knows the baryon fraction that forms 
the central black hole of each quasar and some other astrophysical parameters. Haehnelt 
[59] takes 1% as a reasonable estimate of the fraction, and using reasonable values for 



the other parameters estimates that <7q(M = 1O 13 M ) > 1.30. Multiplying the baryon 
fraction by 10 (or making equivalent changes in the other parameters) reduces this to 
<7q{M = 1O 13 M ) > 1.25, which is shown in figures 2 and 3.0 One sees that this is not very 
constraining. He also gives bounds for lower masses, down to M = lO n M but they are 
even less restrictive. (A more recent estimate (5(| using newer data suggests that a better 
constraint might now be possible, but a detailed analysis has not yet been done.) 

More restrictive bounds fl8| come from damped Lyman alpha systems. The data 
indicate |30| that the baryons in such systems account for f£i ya i p h(> M, (3.0 ± 0.5)) = 
(0.0029 ± 0.0006)/i~ 1 . Let us take the lower bound of the redshift and O, to obtain a lower 
bound on cjq (we set h = 0.5 as its uncertainty is not significant). Then if the baryon frac- 
tion is the same as the universal fraction ~ 0.05, we deduce that Q(> Mi ya i p h, 3.0) > 0.092 
and (Jo(R) > 4.0 where R is the scale corresponding to Mi ya i p h. To obtain a reasonably firm 
bound let us divide this result by 10 (corresponding, say, to assuming that the baryons are 
concentrated by a factor 10). This leads to ao(R) > 2.6. 

The corresponding bound on o"o(i?)/o"o anon (i?) depends on Mi ya i p h, but the dependence 
is much weaker than it is for clusters because 0Q anon (.R) is varying more slowly with scale. 
The bounds displayed in Figures 2 and 3 are for Mi ya i p h = 3 x 1O 11 M , corresponding to 



16 For the purpose of this article 'quasar mass' is taken to mean the total mass of the object concerned 
though the term 'quasar' usually refers only to the region around the central black hole. 

17 The weak dependence on the baryon fraction arises because fiquas(> 1O 13 M0,3) is tiny, so that the 
prediction for it is an exponentially decreasing function of ao 



R = 0.50/1" 1 Mpc. Going down to 1O 1O M only lowers a (R) / a^ non (R) by 17% and going 
up to 10 12 Mq only raises it by 7%, the corresponding scales being R = O.Wh^ 1 Mpc and 
R = 0.76/i -1 Mpc. None of these changes is very significant. 

Our Lyman alpha constraint corresponds to a constraint £l v < 0.2 (with n = 1 and 
h = 0.5). This is the same as that obtained in most (but not all) previous investigations 
[18, p2], which have used both the Press-Schechter formula and numerical simulations. 



6.4 Other observations 

The observations that we have considered are the only ones that can be related more or 
less directly to the linearly evolved density contrast. Other observations can in principle 
yield information after comparing them with numerical simulations, but at present the 
uncertainties seem to be too big to make the additional information useful. For example, the 
pairwise galaxy velocity dispersion has been widely compared with numerical simulations 
and the received wisdom is that with the COBE normalization the prediction is significantly 
too high (eg. ju]]). But others |3l| argue that when the observational (6^] and theoretical 
uncertainties are taken into account the discrepancy ceases to be significant. 



6.5 Summary of the data on a (R) 

The most accurate data point is the one coming from the COBE measured CMB anisotropy, 
which is indicated schematically by the point at 4000/i~ 1 Mpc. That is, of course, the reason 
why we have used it to normalize the theoretical curves. 

The other points all come from galaxy and cluster observations, and for clarity we will 
focus on <7o(-R)/cro anon (-R) as displayed in Figure 3. Most of the points come from the number 
density contrasts of various types of object, as analysed by Peacock and Dodds, which are 
represented by open squares. Within the box each point can be moved down more or less 
independently, and the whole box can be moved up about 20% as indicated by the outer 
error bars. (Both of these uncertainties are more or less at the 1-a level, though the 20% 
estimate may be somewhat high according to this criterion.) The most striking feature of 
these points is the slope of cr (R) / 'crg anon (R) . In the regime Wh' 1 Mpc < R < 40/j" 1 Mpc it 
is positive and rather well defined. This positive slope is seen in both the galaxy correlation 
function (in many different surveys) and in the cluster correlation function [63], and in 



addition it seems to be needed to make N-body simulations give the correct normalization 
for the cluster correlation function 62]. In other words, it is rather firmly established. On 



smaller scales the slope is flatter, and there is also a hint of flattening on larger scales. 

What about the overall normalization of these points? The 20% uncertainty indicated 
by the outer error bars corresponds to at at least a 1-a confidence level in the context of the 
Peacock and Dodds analysis, and both the upper and lower limits do have confirmation from 
other types of data. The lower limit corresponds to the upper limit of Peacock and Dodds' 
bias factor determination bj = 1.0±0.2, which is confirmed by the result [54] bj = 0.7^02 a t 
95% confidence level, that comes from combining the galaxy distribution with observations 
of the bulk flow. The lower limit of the normalization is also confirmed, less strongly, by the 
bulk flow data on the scale R ~ 90/i -1 Mpc. Coming to the upper limit on the normalisation 
of cro(i?)/o"o anon (i?), one sees from Figure 3 that it is confirmed on the scale R ~ 10/i -1 Mpc 
by the abundance of present day rich galaxy clusters if one assumes that their mass is less 
than 2 x 1O 15 M0. A mass in this range is certainly commonly assigned to these rich clusters, 
and indeed a lower (median) mass corresponding to a lower upper limit on o"o(8/i _1 Mpc) 
has been advocated |57p . Essentially the same upper limit is confirmed by the abundance 
of z = 0.3 clusters as cited earlier. 

All this is on scales 10 to 100 Mpc. A crucial lower bound on the normalization on 
some scale in the decade 0.1 to IMpc is provided by the abundance of damped Lyman 



alpha systems at z = 3. In contrast with the case of galaxy clusters, the uncertainty about 
the precise scale matters little here, because the theoretical curves that one is trying to 
normalise are almost flat. It is noteworthy that this lower limit is already presaged by the 
flattening of the slope seen in the regime 4/i _1 < R < Wh- 1 Mpc. 



6.6 Theory versus observation 

The first conclusion to draw from the comparison of theory and observation is that the 
canonical CDM model does not work. As discussed in the Introduction three fixes are 
possible. One is an 'old age' model where h is very low. Another is a 'tilted' model where 
n is significantly below 1. The third is to alter the nature or the amount of the cold dark 
matter, and the specific option that we have explored is to introduce some fraction £l u of 
hot (neutrino) dark matter while keeping the total matter density fixed at the critical value. 

Figure 3 shows that any of these three options is capable of giving an acceptable fit to 
the data, the required parameters being h ~ 0.3 for the first option, n ~ 0.7 for the second 
and f2„ ~ 0.15 to 0.20 for the third. In each case more extreme values (smaller h or n or 
bigger VL V ) are ruled out be the observational lower bounds on the normalization of o~o(R), 
especially the one coming from damped Lyman alpha systems (recall that the one we show 
is supposed to be rather firm). Less extreme values are ruled out by the slope of cro(R) (and 
in the case of the MDM option also by the upper limit on its magnitude) in the tens of Mpc 
range. 

What can we say about the relative merits of the three options? As noted earlier, in 
most (particle physics motivated) models of inflation that implement the tilted option there 
is also a gravitational wave contribution whose effect is to lower the COBE normalization by 
[1 + 6(1 — n)] -1 / 2 . For n = 0.7 this factor is 0.6 which leads to a gross contradiction with the 
data. As a result the tilted option is not viable in most models of inflation, the only known 



exception being 'natural' inflation [64, 15]. The old age option is also problematical because 
of the low Hubble constant that it invokes, which leaves the MDM option as perhaps the 
most attractive of the three. 

So far everything is quite simple, but we have only allowed the three parameters to 
depart from their canonical values h = 0.5, n = 1.0 and £l u = one at a time. We would 
like to understand the entire parameter space. In particular we would like to know the 
observational constraint on n because, especially in conjunction with gravitational waves, 
it is such a useful discriminator between models of inflation. 

One way forward |2(], |2l], ^3| is to perform a least squares fit, though this is somewhat 
problematical because the damped Lyman alpha systems provide only a bound, and because 
it is not clear that all of the points coming from galaxy and cluster correlations are inde- 
pendent. An alternative Q is to simply demand that the curves go through selected points, 
and have slopes within the limits indicated by the box. These options will be explored in 
a future paper pq| , but for the moment we point out that rather definite conclusions flio|| 
about n can already be drawn just from Figure 3. 

To start with, note that in the regime where we have data, a change in h is roughly 
equivalent to a change in n according to the formula Ah = —An. This can be seen by 
comparing the h = 0.3 and n = 0.7 curves in Figure 3, especially for the crucial slope. 
Accepting this and taking 0.4 < h < 0.6, one concludes that with critical density CDM 
requires n ~ 0.6 to 0.8. 

The crucial point now is that all known modifications of the pure, critical density CDM 
hypothesis reduce the COBE normalized o~o(R), because CDM maximises the growth of the 
density perturbation on small scales. Thus the lower limit n > 0.6 will hold in any version 
of the standard model. If gravitational waves reduce the COBE normalization by a factor 
[1 + 6(1 — n)] -1 / 2 , it is not difficult to verify that the lower limit is tightened to n > 0.8 

What about the upper limit? On scales R < 40/i _1 Mpc or so one can clearly cancel an 



Table 1: Estimates of the re-ionisation epoch 1 + z\. 
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increase in n by increasing the fraction of hot dark matter, but at least within the critical 
density MDM option this is impossible on larger scales. Assuming that no cancellation is 
possible, the upper bound n < 1 indicated by the four largest scale galaxy /cluster correlation 
points shown in Figure 3 can only be evaded by decreasing h, which (assuming that h > 0.4) 
allows it to go up to 1.1. A similar upper bound comes from the fact that too much hot dark 
matter gives a strong violation of the damped Lyman alpha system bound, which cannot 
be cancelled by increasing n because of the data in the tens of Mpc range. The conclusion 
is that, at least within the MDM option, there is an upper bound n < 1.1 or so. 

The epoch of re-ionisation 

We end with a brief discussion of estimates of the re-ionisation epoch Z[ on , which are as yet 
rather preliminary. 

Observationally, a lower limit of order z- lon > 5 is well established [56|, and a useful 



upper limit or even a value may soon become available from the small scale cmb anisotropy 



To estimate the re-ionisation epoch theoretically, one starts by estimating the fraction 
/ of the mass that has to collapse into galaxies before re-ionisation occurs. At present the 
astrophysics is not understood sufficiently well to allow a firm calculation of /. Tegmark, 
Silk and Blanchard [S7| give a 'best guess' / ~ 4 x 10 -3 with upper and lower limits 
2 x 10~ 5 < / < 0.4. 

If M m j n is the mass of the smallest galaxies existing at the epoch of re-ionisation, then 
Q(> M mm , Zi on ) = f, and the Press-Schechter estimate gives z- lon in terms of 0"o(-R mm ) and 
/. We take as a best guess M m j n = 10 6 M Q (the same as the minimum mass of present 
day galaxies), while noting that theoretical estimates vary from 10 5 M© to 1O 8 M [§?]]. The 
corresponding scale is R m i n — 0.01 Mpc, two orders of magnitude smaller than the scales 
that we have considered so far. 

If / ~ 1 re-ionisation will not occur until the epoch of non-linearity given by (1 + z Q \) = 
0"o (-Rmin)) but if / is smaller it will be earlier. Tegmark et a/'s 'best guess' gives (1 + z\ on ) ~ 
1.7(7n (-R min ) whereas the lower end of their allowed range gives (1 + z- lon ) ~ 2.5<7o(-R mm ). 
The range is not too big, because if / is small it is an exponentially decreasing function of 

00 (-Rmin)- 

Tegmark et al go on to calculate o"o(0.01 Mpc)/o"o(8/i _1 ) in various models. Inserting 
the values of ao(8h~ 1 ) following from our COBE normalization we deduce the results in 
Table 1. The low estimate of z- lon takes it to be the epoch of non-linearity, corresponding 
to / ~ 1, the other two correspond to the values of / mentioned earlier. In the third row 
we have allowed for the slower growth of the MDM density contrast in the same way as 
Tegmark et al (following |l{J), and in the last row we have used the scaling cfq{R) oc h 2 that 
is appropriate in the small scale regime where cro(R) is practically independent of R. Note 
that the third row is for 30% hot dark matter, whereas at most 20% or so is viable. Thus 
a viable MDM option will give earlier re-ionisation than the other two options. 

Although the uncertainties are big, it seems clear that viable versions of the standard 



model do not lead to very early re-ionisation. On the other hand a re-ionisation redshift in 
the low tens, which might be big enough to be affect the cmb anisotropy, is not out of the 
question. 



7 Conclusion 

To conclude, we have provided a review of the machinery required in order to translate an 
initial spectrum of density irregularities into a form amenable for fairly direct comparison 
with a range of observations. Although some of our discussion has wider application, we have 
focussed on the standard model of structure formation, which invokes an initially Gaussian, 
adiabatic, scale invariant density perturbation and more or less cold dark matter. It is at 
present the one under the most active consideration, because it is the simplest, and also 
because the required initial density perturbation might be generated as a vacuum fluctuation 
during inflation. 

We have gone on to compare the standard model with presently available data, indicating 
the relative merits of altering the Hubble constant, tilting the primordial spectrum and 
incorporating a component of hot dark matter. All of these have been envisaged as ways 
to remedy the shortcomings of the canonical cold dark matter model. 

Our treatment of the observations breaks new ground in that it combines for the first 
time several relatively new observational constraints. For the normalization of the spectrum 
from the COBE measurement of the cmb anisotropy, we use the most recent determina- 
tion that has been provided by Gorski and collaborators. It is significantly higher than 
earlier determinations. We assess carefully the uncertainty in the determination of the nor- 
malization on the scale 8/i _1 Mpc that comes from the cluster abundance, quantifying its 
dependence on the assumed cluster mass. We establish a simple and rather firm constraint 
from damped Lyman alpha systems at high redshift, which provide a powerful lower limit 
on some scale of order 0.1 to IMpc. Finally, on the basis of an earlier study [67] we make 
a preliminary estimate of the epoch of re-ionisation in the various models. 

Our technique for comparing theory with observation is to look not at the power spec- 
trum, but at the dispersion of the density contrast smoothed on a scale R. This has the great 
advantage of being closer to what is actually measured, while still being simple to calculate 
for a given theoretical model. Particularly when one normalises theories and observations 
to a benchmark model as in Figure 3, the comparison between theory and observation can 
be very clearly illustrated. 

While present observations are strong enough to exclude the standard CDM model, they 
are not of sufficient quality to select amongst different ways of generalising it, particularly 
when one realises that if tilt or a hot component are to be varied then one must certainly 
also allow h to vary in combination. It is clear from Figure 3 that these options are most 
cleanly probed by observations on scales from a few tens to a few hundreds of megaparsecs, 
and we should look forward to this region of the spectrum being probed by both small 
scale microwave anisotropy experiments and by larger scale galaxy correlation (and peculiar 
velocity) measurements. 
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Figure Captions 

Figure 1. [NOT AVAILABLE VIA BULLETIN BOARD] 

This figure is reproduced from []5C| ] , and it shows predictions for the spectrum Vt(1) of the 
cmb anisotropy. The full line 'S' shows the contribution (normalized to 1 at I = 2) of an 
adiabatic density perturbation with n = 0.85, h = 0.5, Qb = 0.05 and critical total density. 
It was calculated for pure CDM, but it is insensitive to the nature of the dark matter. (The 
dashed line is the prediction with Qb = 0.01 and the same value for h, but note that these 
quantities are actually linked by the nucleosynthesis relation O^/i 2 = 0.013±.002.) The line 
'T' shows the gravitational wave contribution, also normalized to 1 at I = 2. As discussed 
in the text, it is not yet known whether this contribution is actually significant. 

Figures 2a and 2b. 

Observational data, shown in comparison to the canonical CDM model plus four variants, 
two with a hot component, one tilted and one with a low Hubble constant. All of the models 
are normalized on large scales to the COBE data as described in the text, and this data is 
represented schematically by the filled triangle at Rh -1 ~ 4000 Mpc. The rest of the data 
comprises damped Lyman alpha systems and quasars (two lower limits on the left hand 
side), cluster abundance (triangle, with a tilted error bar indicating the dependence of the 
prediction on the assumed cluster mass), galaxy correlations (squares, error bars explained 
in text and in figure 3 caption) and bulk velocities (star). The second figure is a closeup of 
the first, which omits the COBE point. 



Figure 3. 

This is the same as Figures 2a and 2b except that everything has been divided by the 
prediction of the canonical CDM model. As explained in the text the galaxy correlation 
points enclosed by the 'error box' can be moved up and down more or less independently 
within the box, and in addition the whole box can be moved up and down by the amount 
indicated by the outer error bars. With this exception all of the error bars are more or less 
statistically independent, and represent something like a 1-a confidence level. The lower 
bounds are intended to express a high level of confidence as explained in the text. 
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